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Abstract. We study in the BFKL approach the total hadronic cross section for the collision of 
two virtual photons for energies in the range of LEP2 and of future linear colliders. The BFKL 
resummation is done at the next-to-leading order in the BFKL Green's function; photon impact 
factors are taken instead at the leading order, but with the inclusion of the subleading terms required 
by invariance under changes of the renormalization scale and of the BFKL scale sq. We compare 
our results with previous estimations based on a similar kind of approximation. 
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INTRODUCTION 

The total hadronic cross section for the collision of two off-shell photons with large 
virtualities is a fundamental observable, since it is fully under the control of perturbative 
QCD. It is widely believed that this total cross section is the best place for the possible 
manifestation of the BFKL dynamics y]] at the energies of future linear colliders (for a 
review, see Ref. [2]]). For this reason, many papers [3] have considered the inclusion 
of the BFKL resummation of leading energy logarithms. In a remarkable paper [0] 
(see also Ref. ^\), BFKL resummation effects have been taken into account also at 
the subleading order and evidence has been presented that the appearance of BFKL 
dynamics is compatible with experimental data already at the energies of LEP2 [@,0]. In 
this work [8] we estimate the energy dependence of the 7*7* total hadronic cross section 
in an energy range which covers LEP2 and future linear colliders. The procedure we 
follow is approximate, since we use the singlet forward NLA BFKL Green's function 
together with forward 7* — > 7* impact factors at the leading order. However, in the 
impact factors we include the subleading terms required by the invariance of the full 
amplitude at the NLA under change of the renormalization scale and of the energy scale 
so entering the BFKL approach. 



THE 7*7* TOTAL CROSS SECTION: NUMERICAL ANALYSIS 

The total hadronic cross section of two unpolarized photons with virtualities Q\ and 
Q2 can be obtained from the imaginary part of the forward amplitude. Following the 
procedure of Refs. |@], it is possible to write down the cross section with the inclusion 



of NLO corrections in the Green's function only, while keeping the impact factors 
at the LO. In fact, the requirement of invariance of the amplitude at the NLA under 
renormalization group transformation and under change of the energy scale sq allows to 
fix the Hr- and so -dependent terms in the NLO impact factors 

A(*)=z(v)m^, «w=t ta i- (1) 

The details of the analytical calculation can be found in Ref. BSfl. We use the series 
representation [9], that is one of the infinitely many possible ways, equivalent with NLA 
accuracy, to represent the total cross section. In the case of the y*y* — > VV process 
where V stands for light vector meson (p, CO, 0), it turned out that the contribution to 
the amplitude from the kernel starts to dominate over that from the impact factors in 
the series from n = 4. This makes evident the fact that the high-energy behavior of the 
amplitude is weakly affected by the NLO corrections to the impact factor. Therefore, 
our approximated y*y* total cross section should compare better and better with the 
correct result as the energy increases. In order to stabilize the perturbative series, it is 
necessary to resort to some optimization procedure, exploiting the freedom to vary the 
energy parameters, /i^ and sq, without corrupting the calculation but at the next-to-NLA. 
We use both the principle of minimal sensitivity (PMS method) [10] and the Brodsky- 
Lepage-Mackenzie (BLM) method 111 ill : for some selected values of the energy s in 
the region of interest the optimal scales ji R and sq are found and the cross section is 
thus determined. Then, the curve giving the cross section vs the energy is obtained by 
interpolation. 

In order to compare the theoretical prediction with the existing data from LEP2, 
we cannot neglect the contribution from LO quark box diagrams lfl2ll which is of 
order a 2 (lns)/s. On the other hand, the soft Pomeron contribution, if estimated within 
the vector-dominance model, is proportional to Cy*y* ~ (my/Q 2 ) 4 Oyy and is therefore 
suppressed for highly virtual photons. We restrict ourselves to the case of a symmetric 
kinematics, which means the same virtuality <2i = Qi = Q for the two photons. This 
is the so-called "pure BFKL regime", as opposite to the "DGLAP regime" realized for 
strongly ordered photon virtualities. In Fig. [TJ (left) we summarize our results for the 
CERN LEP2 region: we show the NLO BFKL curves obtained by the PMS and the 
BLM methods, to which we added the contribution of the LO quark box diagrams. 
For comparison we put in this plot also the experimental data from CERN LEP2, 
namely three data points from OPAL [|7|] (<2 2 =18 GeV 2 ) and four data points from 
L3 [@] (<2 2 =16 GeV 2 ). We observe first of all that the difference between the two 
theoretical curves can be taken as an estimate of the systematics effects which underlay 
the optimization procedures adopted here. Then, the comparison with experimental data 
is acceptable, although within uncertainties which are large both on the theory and on the 
experiment side. Around the energy for which the condition for the BFKL resummation, 
a s Y ~ 1, is satisfied, which in the present conditions corresponds to Y ~ 5, both the 
PMS and BLM curves agree with experimental data within the errors. Finally, we remark 
that the determination from Ref. JUl falls between our two curves from PMS and BLM 
methods. From the energy dependence of the NLO BFKL cross section determined 
through the PMS method at <2 2 =17 GeV 2 we obtained also the "dynamical" Pomeron 
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FIGURE 1. (Left) Energy dependence of the total cross on section on Y = \n(s/Q 2 ). For comparison, 
experimental data from OPAL |7] (stars, Q 2 =18 GeV 2 ) and L3 [(fl (diamonds, Q 2 =16 GeV 2 ) are shown. 
(Right) F-dependence of the Pomeron intercept (minus 1) calculated from the total cross section with the 
PMS method. All theoretical curves are obtained for Q 2 =17 GeV 2 and rif = 4. 



intercept (minus 1) as a function of the energy. 

The result is shown in Fig. Q] (right). In Fig. 2 we show the /"-behavior of the total 
cross section for Q 2 =20 GeV 2 (n/=5) in an energy region not explored by past and 
present experiments, but relevant for future colliders. We plot here the two curves 
obtained in the present work with the PMS and the BLM methods. The condition for 
the BFKL resummation, d s Y ~ 1, corresponds here to Y ~ 6; around this energy the 
deviation between the PMS and the BLM methods is about 50%. This discrepancy can 
be taken as an estimate of the systematic uncertainty of this approach. We observe that 
our determination from the BLM method is in quite good agreement with the result of 
Ref. |0] (see Fig. 4 of that paper), obtained for the same kinematics. 



CONCLUSIONS 

We have found that, if suitable methods are used to stabilize the perturbative series, a 
smooth curve for the energy behavior of the cross section can be achieved. Our result in 
the CEPvN LEP2 region compares favorably with experimental data. Systematic effects 
coming from the optimization procedure are estimated by the comparison with two 
different methods. Our findings in the CERN LEP2 region are in agreement with the 
result of Ref. [4]], where for the first time subleading BFKL effects were considered in 
the 7*7* total hadronic cross section. In comparison with the latter work, which deals 
with the same process under consideration in the present paper, the elements of novelty 
are the following: 

• the optimization procedures to stabilize the perturbative series are performed on 
the amplitude itself and not on the NLO Pomeron intercept, which is not a physical 
quantity; 

• the impact factors, although taken at the LO, contain the appropriate NLO terms, so 
that the dependence on the energy scales entering the process (the renormalization 
scale fiR and the parameter sq introduced in the BFKL approach) is pushed to the 
next-to-NLA; 
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FIGURE 2. The ^-dependence of <jy>y> [nb] for virtual photon - photon collisions at (2 2 =20 GeV 2 and 
tif = 5, obtained with the PMS and the BLM methods. 



• two optimization methods are used, thus having a control of systematic effects at 
work. 

The numerical effect of the neglected subleading corrections to the impact factors 
cannot be quantified. We expected that it be modest in the second region of energy 
considered in this work. Here our prediction should be very close to the true NLO result. 

The final word will be said when the y*y* cross section will be calculated fully in the 
next-to-leading approximation. 
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